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Abstract 

Transport current distribution in clean 2i2"-NbSe2 crystals is studied by mea- 
suring the self-induced magnetic field across the sample. Below T c most of the 
current flows at the edges of the crystals due to strong surface barriers, which 
are found to dominate the transport properties and the resistive transition. 
The measured critical current is determined by the critical current for vortex 

penetration through the surface barrier rather than by bulk pinning. 
PACS numbers: 74.60.Ec, 74.60.Ge, 74.60.Jg 
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Vortices in type-II superconductors have to overcome surface and geometrical barriers 
(SB) in order to exit or penetrate into the superconductor. The effect of SB on the mag- 
netic properties has been extensively studied theoretically |I]~§J- in recent years numerous 
studies have shown that SB dominate the magnetization behavior in clean crystals of high 
temperature superconductor (HTSC), in particular at elevated temperatures p|-|TT]|] . Since 
bulk pinning is strongly reduced by thermal fluctuations, the relative importance of the 
SB is expected to grow with temperature 0]. At low temperatures, in contrast, and in 
low-T c superconductors in general, bulk pinning is expected to be the main source of hys- 
teretic magnetization. Recent theoretical works [|T|,|13|] have suggested that SB can also 
significantly modify the transport properties of superconductors. Current distribution mea- 
surements in Bi 2 Sr 2 CaCu208 (BSCCO) crystals |14l , |T5|| have revealed that over a wide range 
of temperatures and fields, the transport current indeed flows predominantly at the edges of 
the crystals where vortices enter and exit the superconductor. In this paper we report that 
SB dominate the transport behavior also in the low-T c superconductor NbSe2- SB should 
therefore be of significant importance in a wide range of superconductors in affecting the 
transport as well as the magnetic properties. 

High purity 2if-NbSe2 crystals with T C =7.2K were grown as described previously ||TB| . 



Several crystals were cleaved and cut into a rectangular strip shape. The data presented 
here are for one of the crystals with dimensions of 1.46 mm(/) x 0.35 mm(w) x 0.04 mm((l). 
The features reported below were observed in all the investigated crystals. Four pads for 
electrical contacts of 0.1 x 0.1 mm 2 with 0.13 mm separation between pairs were prepared 
by Ag/Au evaporation. The crystal was mounted onto an array of 19 two-dimensional 
electron gas GaAs/AlGaAs Hall sensors 10 x 10 /mi 2 each with a 10 /mi separation. The 
inset to Fig. 1 shows a schematic top view of the experimental setup. A dc magnetic field 
Hjc was applied perpendicular to the plane of the sensors and parallel to the c-axis of the 
crystal. An ac transport current I ac in the range of 1 to 30 mA at 65 Hz was applied, and 
the corresponding self-induced magnetic field B ac (x) across the crystal was measured by the 
Hall sensor array. The four-probe resistance of the sample was measured under the same 



conditions using a lock-in amplifier. 



There are three main regimes for the flow of the transport current as described in Ref. |I3 
(a) Uniform current flow in the bulk of the sample. In this case, applying the Biot-Savart 
law, the perpendicular component of the self-induced field, as measured by the sensors, 
decreases monotonically from the left edge to the right, (b) Surface barrier dominated flow 
for which most of the current flows at the two edges of the crystal in order to drive the 
vortices over the entry and exit SB. The resulting self- induced field across the crystal is 
opposite in sign and increases from left to right, (c) Vortex pinning which prevents vortex 
motion and results in zero self- field inside the sample. The corresponding transport current, 
in this case, has a characteristic Meissner distribution with a reduced current in the center 
and enhanced at the edges |L7]-[n|. 

Figure la shows an example of the self-induced field B ac as a function of temperature at 
ifd c =0.1T and J ac =6mA. For clarity, only the sensors 3 to 18, which are under the crystal, 
are shown. The corresponding resistance is shown in Fig. lb. The observed behavior is very 
similar to that reported for BSCCO crystals |H |. Above T c , the current flows uniformly 



as expected, and B ac decreases monotonically from the left edge of the crystal (sensor 3) 
to the right edge (sensor 18). This behavior is seen more clearly in Fig. 2a, which shows 
B ac {x) profile for all the sensors at 7.6K. As the temperature is decreased, the SB sets in 
immediately below T c drawing an increasingly larger fraction of the current to the edges. 
This gives rise to the observed temperature dependence of B ac in Fig. la: a drop in B ac 
just below T c , followed by crossing of the curves and sign reversal at lower temperatures. 
At the crossing point, half of the transport current flows across the bulk and the other half 
at the two edges of the sample. At T = T max (see Fig. la) most of the current flows at 
the crystal edges and the self-field profile inside the crystal is completely inverted relative to 
the uniform flow shown in Fig. 2a for T=6.8K. Note that any bulk mechanism or 

finite skin depth effect will result in B ac which is either zero (perfect shielding) or positive 
in the left half of the sample (finite shielding), but it cannot cause B ac to become negative 
(negative permeability [[|I|). The sign reversal of B ac is a unique property of the SB. As 



the temperature is further decreased, the vortices become immobile for T < Td due to the 
combination of the SB and bulk pinning, resulting in a vanishing B ac response within the 
crystal as shown by the 4.9 K profile in Fig. 2a. 

In Ref. [|nj arrays of seven Hall sensors were used, which allowed mapping of B ac typically 
over only half of the sample width. Here we have extended the arrays to 19 sensors that 
provide more detailed B ac (x) over the entire sample width. This improvement has two 
significant advantages. First, we can readily examine both edges of the sample and evaluate 
the symmetry of the current distribution. The second major advantage is that having 
19 values of B ac field across the sample is sufficient in order to directly invert the field 
distribution into the current distribution. We represent the sample by 19 current filaments 
located at about half of the crystal thickness and equally spaced across the width. The 
currents in the filaments are obtained by inverting the 19 x 19 matrix that transforms 
between the current and the field values using the Biot-Savart law. Figure 2b shows the 
obtained I(x) corresponding to the three field profiles in Fig. 2a. As expected, above T c the 
current flows uniformly across the crystal (7.6K profile). However, below T c the current starts 
to accumulate at the edges due to the strong SB. At T = T max practically all the current 
flows at the two edges in a form of two 5 functions with negligible current in the bulk (6.8K 
profile). At low temperatures, 4.9K profile, the vortices become immobile and the current 



is distributed in the corresponding self-shielding form [p^ -[T9| with a minimum in the center 
and a rapid increase near the edges. Note that the actual current distribution is continuous 
across the width and the thickness of the crystal, and hence our derivation of the discreet 
current filaments in Fig. 2b is only an approximation. Yet, this simple analysis clearly 
visualizes the underlying mechanisms, and in particular the main transition from a uniform 
current flow above T c to SB dominated flow below T c . This finding shows that in clean 
NbSe2 crystals the vortex flow rate is determined by the transmission probability through 
the SB rather than by bulk vortex dynamics. An important experimental implication is that 
the transport measurements in this case reflect the resistive properties of the SB rather than 
the bulk properties. 



Another aspect of the SB is its asymmetry with respect to vortex entry and exit [|TJ. 
Vortex entry requires a larger force than vortex exit. As a result a larger current flows at 
the vortex entry edge in order to maintain the same vortex flow rate throughout the sample. 
The role of the edges is interchanged as the direction of the ac current changes, with a larger 
current flowing on the opposite edge of the crystal. This mechanism results in a significant 
local second harmonic self-field signal, as shown in Fig. 3 for i^ c =0.5T and J ac =6mA. 
Second harmonic is a unique feature of the SB due to its asymmetry with respect to the 
current direction. Bulk vortex dynamics, in contrast, results only in odd harmonics, since 
bulk I-V characteristics are antisymmetric with respect to the current. Figure 3 shows that 
the second harmonic signal, and hence the SB, set-in immediately below T c concurrently 
with the resistive drop. The narrow dip in the resistance at T p in Fig. 3 is the common peak 
effect in NbSe2 p0|-f23|, which in our high purity crystals is extremely narrow (see also Fig. 
la). Within our experimental resolution we do not observe substantial changes in current 
distribution in this narrow peak-effect region. 

In contrast to BSCCO, the thermal activation of vortices in NbSe2 is weak due to much 
lower T c and lower anisotropy. We can therefore simplify the description by analyzing 
the behavior in Fig. 3 in terms of bulk and surface barrier critical currents, l\ and J*, 
respectively. Let us first consider the case of finite bulk pinning with no SB, namely /,?=0. 
The bulk critical current l\ is generally expected to increase with decreasing temperature, 
but as long as the applied current is larger than l h c iT\ the current should flow uniformly 
across the sample [|17HT9|1. As a result, B ac should be positive like above T c , and there should 
be no second harmonic signal. When I^(T) approaches the value of the applied current at 
some characteristic temperature T max , vortices stop moving, and B ac should drop rapidly 
from the full positive value to zero as the temperature is further decreased. This scenario is 
inconsistent with the data in Fig. 3a: B ac is negative instead of positive above T max , and in 
addition, significant second harmonic is present. Furthermore, below T max , B ac does drop 
rapidly, but this drop occurs while B ac is negative rather than positive. 

We now consider the opposite scenario of finite SB with no bulk pinning, 1^=0. Since 
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I* CO is related to the critical field H C (T), it is expected to grow rather linearly below T c . 
Therefore, upon cooling, a progressively larger fraction of the applied current should be 
drawn towards the edges. The corresponding B ac should drop approximately linearly from 
a full positive value at T c down to a full negative value at a characteristic temperature at 
which J* (T) reaches I ac , and as a result practically all the applied current flows at the edges. 
In this temperature interval a significant second harmonic signal should be also present, 
as explained above. The corresponding resistivity of the sample should decrease gradually 
below T c with a sharp drop towards zero at the same characteristic temperature. The data 
in Fig. 3 above T max is fully consistent with this scenario, thus allowing us to identify T max 
as the characteristic temperature at which I ac = I^{T). 

Now let us examine the behavior below T max where I ac < /f (T) and no vortices can exit 
or penetrate into the sample (in absence of thermal activation). Yet, if 1^=0 the vortices 
can move freely inside the sample and change their distribution according to B ac imposed 
by the full I ac flowing on the edges. Since the edges have already absorbed all of I ac , there 
should be no further change in the current distribution below T max , and B ac should remain 
constant at its full negative value (as observed in BSCCO jnj). O n the other hand, if a 
small bulk I h c is present in addition to J*, B ac should drop sharply to zero below T max since 
vortices become immobile. T max in this case corresponds to the temperature at which I ac 
becomes equal to the total critical current I C {T) = I* + 1^. The pronounced drop of B ac at 
T max in Fig. 3a indicates therefore the existence of some finite bulk I h c . The fact that B ac is 
practically fully inverted at T max shows, however, that l\ <C l s c . Note also that the second 
harmonic signal should be absent below T max , since no vortices penetrate through the SB, 
consistent with the data in Fig. 3a. The tail of B ac as well as the weak resistive tail below 
T max indicate, furthermore, that a small flux creep is present resulting in some finite vortex 
motion. Note, that even in this weak creep regime B ac is inverted, showing that most of the 
current flows at the edges and vortex dynamics is governed by SB. The vortices become fully 
immobile only at a lower temperature Td, below which the creep stops within our resolution. 

Finally, we demonstrate here that the above condition, that the total critical current I c is 
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determined mainly by the critical current of the SB, holds over the entire investigated range 
of temperatures and fields. At low temperatures and relatively low current, the vortices 
are immobile and the current flows with the characteristic Meissner distribution resulting 
in vanishing B ac . In this case one concludes that the applied current I ac is below the total 
critical current I c — l\ + I*, but the values of the individual critical currents cannot be 
determined. In order to gain this important information one has to increase the applied 
current just slightly above the total critical current. In this situation the applied current is 
precisely divided between the bulk and the edges according to l\ and I s c . As a result the 
vortices are set in motion, and the corresponding B ac signal provides the information on the 
relative importance of l\ and l s c . Such a measurement is presented in Fig. 4 which shows 
B ac , as measured by sensor 3, at various I ac between 5 and 30 mA at 0.3T. The arrows 
indicate the temperatures T max at which I ac = I C (T) — l h c + I* for the various applied 
currents as described above. The corresponding I C (T) dependence is shown in the inset. We 
find that in all cases when I ac reaches I C {T) the corresponding B ac signal is fully inverted (as 
the 6.8K profile in Fig. 2a) indicating that practically all the current flows at the edges, and 
that l\ <C Thus, in clean crystals of NbSe 2 the measured I c reflects mainly the critical 
current of the SB and not the bulk critical current. 

In summary, the distribution of transport current in clean platelet crystals of 2if-NbSe2 
was studied by measuring the self-induced magnetic field. The use of extended arrays of 
Hall sensors allows for a direct inversion of the self-field profile into the current distribution 
profile across the crystal width. Below T c the current is found to flow predominantly at the 
sample edges due to strong surface barriers. The SB govern the apparent resistivity of the 
crystals. Furthermore, the measured critical current is determined by the critical current of 
the SB barrier rather than by bulk critical current. 
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FIGURE CAPTIONS 

Fig. 1. (a) Transport current self-induced field B ac as measured by sensors 3 to 18 across 
the width of a NbSe2 crystal at Hd c = 0.1T and I ac = 6mA. At the crossing point of the 
curves half of the transport current flows in the bulk and half at the sample edges. At 
lower temperatures the inverted curves indicate that most of the current flows at the edges. 
Vortices are immobile below T d . (b) Corresponding four probe resistance measurement. The 
resistance drops sharply at T max and shows a tail at lower temperature due to weak flux 
creep. Inset: schematic cross section of the sample mounted on an array of 19 Hall sensors. 

Fig. 2. (a) Three B ac {x) profiles at different temperatures taken from Fig. la. (b) 
Transport current distributions obtained by inversion of B ac in (a). Above T c the current 
flows uniformly across the crystal (7.6K profile). At T —T max — 6.8K practically all the 
current flows at the two edges of the sample due to strong SB. At low temperatures the 
vortices are immobile and the current is distributed in the corresponding self-shielding form 
(4.9K profile). 

Fig. 3. First and second (multiplied by 3) harmonics of B ac as measured by sensor 3 (a), 
and the corresponding resistance (b) at H dc = 0.5T and I ac = 6mA. The SB, and hence the 
second harmonic signal, set in immediately below T c . At T max , I ac —l\ + and B ac , the 
second harmonic signal, and the resistance drop sharply. The narrow dip in the resistance 
at T p is the peak effect in NbSe 2 . 

Fig. 4. B ac as measured by sensor 3 for various I ac between 5 and 30 mA at H dc = 0.3 T. 
As the current is increased vortices become mobile at progressively lower temperatures. The 
arrows indicate the temperatures at which the corresponding I ac equals I C (T). The negative 
B ac values at these temperatures show that I h c <C I s c . Inset: temperature dependence of the 
total critical current I c at H dc = 0.3T. 
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Fig. 1 , Paltiel etal. 
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Fig. 3, Paltiel etal. 
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